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Abstract 



The electromagnetic corrections to the low energy scattering amplitude involving 
charged pions only are investigated at leading and next-to-leading orders in the two- 
flavour chiral expansion. As an application, the corresponding variation in the strong 
2S — 2P level shift is evaluated. The relative variation is of the order of 5%. 
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1. The study of hadronic atoms has become a very active field (see ]1J for a recent account). 
Many experiments are devoted to measure the characteristics of such atoms with high preci- 
sion 0, |, [5L R]. These experimental results carry a well founded theoretical interest, since 
they provide a direct access to hadronic scattering lengths, leading, in this way, to valuable 
informations concerning the fundamental properties of QCD at low energy. For instance, the 
presently running DIRAC experiment aims at measuring the pionium lifetime r with 10% accu- 
racy . This would allow one to determine the difference a[] — with 5% precision by means 
of the Deser-type |7J] relation [§] 

r- l <x (a°-al) 2 , (1) 

where a\ denotes the /-wave nir scattering length in the channel with total isospin I. On 
the other hand, chiral perturbation theory (ChPT) predictions for the scattering lengths have 
reached a precision amounting to 2% ||. Before confronting the experimental determination to 
the ChPT prediction, it is necessary to get all sources of corrections to the relation ([!]), valid in 
the absence of isospin breaking, under control. In this connection, bound state calculations were 



performed using different approaches, like potential scattering theory [JTOj, pi]] , 3-D-constraint 
field theory fI2|| , Bethe-Salpeter equation and non-relativistic effective lagrangians |T3|, |Hi|j . 
For a review on the subject and a comparison between the various methods we refer the reader 
to [|T6 1 . Within the framework of non-relativistic effective lagrangians, the correct expression of 



relation ([I]) which include all isospin breaking effects at leading order (LO) and next-to-leading 



order (NLO) was given as |17 



r- 1 = i a 3 (AM 2 ± - AM 2 ~ Afe a 2 ) ^ A 2 (1 + K) . (2) 
9 

In the preceding equation, a = e 2 /(47r) stands for the fine-structure constant, A and K possess 
the following expansions [|17| in powers of the isospin breaking parameter k e [a, (m^ — m u ) 2 ] 



A = --ReA t r + o(,), (3) 

OZ7T 

K = (a + 2a 2 ) 2 -^(\na-l)(2a° + a 2 )+o(K). (4) 



The quantity of interest, 

3 



32vr 



Re AL.' 00 = a o ~ a l + hi{m d - m u ) 2 + h 2 a , (5) 



represents the real part of the 7r + 7r~ — > 7r°7r° scattering amplitude at order k, calculated at 
threshold within ChPT to any chiral order and from which we subtract the singular pieces 
behaving like q~ x and lng whith q being the center-of-mass three-momentum of the charged 
pions. The coefficient h 2 was calculated in |18|] at next-to-leading order in the chiral expansion 
while hi = 0{m u + m d ) 



The DIRAC proposal [fj] also mentioned the possibility to measure the strong 25 — 2P energy 
level shift AE strong of the pionium. How this measurement could be performed in practice has 
been discussed in Ref. |I9| . A simultaneous measurement of r and A^trong would allow to pin 



down a[| and separately, since [20] 



A-E s trong OC 2(2° + . (6) 
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Bound state calculation of the isospin breaking corrections to (|g) were done in ]TTJ| using 
potential scattering theory (the main contribution to the total level shift comes from vacuum 
polarization effects, see e.g. pi| ). Non-relativistic effective lagrangian calculations concerning 
AE strong are not available for pionium, but exist for the pionic hydrogen case [p2 |. One may 
thus reasonably expect an expression similar to Eq. (§) 

A£ strong = ~a a M n ±A'(l + K'), 

o 

* = ReAt;. ;+ ~ + o{k) , (7) 

OZ7T 

involving the tt + tt~ — > 7r + 7r~ scattering amplitude A + ~ ;+ ~ . Electromagnetic corrections to this 
scattering amplitude will be calculated at next-to-leading order in the present work. 



2. The elastic scattering process 

7l+( P+ ) + 7T-(p_) -> 7T + (p' + ) + 7T-(p'_) , (8) 

is studied in terms of the Lorentz invariant Mandelstam variables 

s = (p++P-) 2 , t = (p + -p + ) 2 , u=(p+-p'_) 2 , 

satisfying the on-shell relation s + 1 + u = 4M 2 ± . These variables are related to the center-of- 
mass three-momentum q and scattering angle 9 by 

s = A(M 2 ± + q 2 ) , t = -2g 2 (l-cosfl), u = -2q 2 (I + cos 9) . (9) 

Let A + ~' + ~ and A ++]++ denote the respective scattering amplitudes for the process (|J) and for 
the crossed channel reaction 7r + 7r + — > 7r + 7r + . Then, s <-> w crossing is expressed as 

= (10) 

We shall calculate the scattering amplitude ([H]) at NLO including electromagnetic effects. The 
strong sector chiral lagrangian for two-flavour ChPT was constructed in [23]. Treating isospin 
violation of electromagnetic origin requires the extension of ChPT in order to include virtual 
photons. This can be done by building operators in which photons figure as explicit dynamical 
degrees of freedom. The electromagnetic sector of the chiral lagrangian for two-flavour ChPT 
has been discussed at NLO in [18| and [24]. We shall work in the m u = rrid limit and use 
the lagrangian representation of [I8|. At one- loop accuracy, all of the following chiral orders 
are present: p 2 ,e 2 ,p 4 ,e 2 p 2 ,e A . We do however not consider the 0(e 4 ) contributions, which 
include two-photon exchange box diagrams, and which are expected to be smaller than the 
other contributions at the same order. Using Feynman graph techniques, the amplitude (|TTJ) 
can be represented at NLO by the diagrams^ depicted in Fig. [[]. 

The scattering amplitude ([TIT) is conveniently written in the following s *->■ t symmetric decom- 
position 

f <? - M 2 

A+-'+-(s,t,u) = |_^ + b+-;+-( M , m ) + c + -;+-m,«) 

+ e 2 (—)iFv(s)} 2 \ + {s^t}. (11) 



3 If one uses the so-called cr-model parametrization of the pion fields, the diagrams (e) of Fig. || and (b) of 
Fig. H vanish identically. 
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(d) (e) (f) 

Figure 1: The various topologies of Feynman diagrams contributing to the charged tttt scattering 
amplitude at order one loop, but ignoring 0(e 4 ) effects. The full circles appearing in the Born- 
type diagram (f) are made explicit in Fig. [| 




Figure 2: The electromagnetic vertex function of a charged pion to one-loop order. The full 
square takes into account the contribution from the low-energy constants just as the tree con- 
tribution including the effect of wave function renormalization. Diagrams of order 0(e 3 p) are 
discarded. 



4 



In the preceding expression, B + ;+ collects the unitarity pieces arising from the diagrams of 
type (c) and (d) in Fig. [I], 



K s,t,u) 



1 



+ 
+ 



2F 4 
1 

Jl 
1 

2e 2 
F 2 
l_ 
F 2 



{s-M* ) 2 J 00 (s) 



(u -t)(s- 4M 2 ±) + 2sA n + 4 A 2 



J + _( S ) 



4 12 

(u - 2M 2 ± - 2A n )(u - 2M 2 ± - 2A W - Ae 2 F 2 )J + _{u) 



+ — (u - 2M 2 ± - 2A W ) [2(s - 2M 2 ± )G + _ 7 (s) - (u - 2il^ ± )G + _»] 

t — u 
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t + u 



(12) 



where 



A w = M 2 ± - M 2 . 



The expressions for the various loop functions appearing in (|i2[) can be found in |]18 |. The 
function C +_;+ ~ represents the contributions from tadpoles as well as from the strong and 
electromagnetic |TS| low-energy constants 

s-Ml e 2 



C 



+-;- 
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+ 
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F 2 32vr 2 
e 2 M 2 



1 + In 



M 2 



32tt 2 F 2 
1 

48vr 2 F 4 
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[( s -2M 2 ± ) 2 (r 1 + r 2 



2tt 2 F 2 



(a - 2M\ 
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t + u 
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1 A. 
96tt 2 F 4 M 2 



(-3s 2 + 16sM 2 + 2uM 2 - 23M^ ) 



(13) 



In this last expression, we have dropped all contributions of order 0(e 4 ) and beyond. In 
particular, we have expanded all logarithms of the pion mass ratio, 



-TT U 



+ 



The terms /C + ;+ and /C ++;++ involve combinations of the electromagnetic counterterms ki, 

4Z\ - 40Z- 



/c ++; - 



3 

-|3 + 



9 / 
AZ 

~9~ 



) h - —k 2 - 9A; 3 + 4ZA; 4 + 4(1 + 8Z)k 6 + 2(1 - 8Z)k 8 , (14) 



fci 



248Z 
9 



k 2 + 9A: 3 - 20Zfc 4 + 4(1 + 8Z)A; 6 + 2(1 - 8Z)k 8 . (15) 



The low-energy constants U and ki are related to the respective renormalized running couplings 
li(fj) = l\ and klip) = k\ at the scale /i = M v ±, 



I; 



li 



327T 2 



M 2 d 
h + In — £ 
/i 2 



A'' 



<7; 



327T 2 



M 2 , 
fci + In — £ 
/i 2 
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where j{ and Cj are the renormalization group beta-functions whose values are given in p3 
and fl8f , respectively. The last quantity that remains to be defined in ([□]) is Fy, which denotes 
the electromagnetic pion form factor 

(n + (p' + )\J^0)\7r + (p + )) = - ie (p + +p' + rF^Q 2 ) , Q = p' + - p + , 



with being the electromagnetic current. At the order we are working, it is given by |23|, B5 



1 + 



p 



2 r 



6F 2 



1 - 



4M 2 



J + -{P 2 ) + 



1 



16vr 2 



(16) 



Finally, let us mention that the scattering amplitude (|TTD is scale independent but infrared 
divergent. This infrared divergence is signaled by the presence of a In (m 2 /M 2 ± ) term, where 
m 7 is a small photon mass acting as an infrared regulator. These infrared divergencies have to 
cancel in observable quantities, like the total cross section with soft photon emission. These 
divergencies do however not show up in the discussion of the electromagnetic corrections to the 
strong level shift, to which we now turn. 



3. The scattering lengths a[j and a 2 , are well-defined quantities in the absence of radiative 
corrections. Their NLO expressions were derived in p3 and, for historical reasons [ 26fl , we shall 
reproduce them in terms of the charged pion mass 



7M 2 ± 

7T 

32ttF 2 

7T 

16ttF? 



5M 2 ± 

7T 

84tt 2 F 2 

2 



h + 2l 2 h H h H 

8 10 8 



M 2 



T ^i" ^ 7- 3 3 

- h + 2l 2 h h + - 

12vr 2 F 2 2 8 2 8 



(17) 



with F n being the pion decay constant whose NLO expression is given in the isospin limit by [53 

F, = F ( 1 + k 



16tt 2 F 2 

We are interested in electromagnetic corrections to the strong energy level shift. To this end, 
one has to consider the part of the amplitude corresponding to one-photon-irreducible (OPI) 
diagrams. It is then convenient to subtract the contribution from diagram (f) of Fig. [I] as 
follows 



.4 



A 



Born 



(s,t,u) 
'(s,t,u) 



A 



s,t,u) + A] 



-=f;+=F/ 



u-t r „ w Nl2 u - s [nI ,. l2 



s 

s - t 



[F V (S)Y + 



[Fy(t)Y 



^Born ~ 

Next, we expand the real part of AqpI +t in powers of q, 
Re AqpI + ~(s, t, u) 
[s,t,u) 



[Fy{n)f+ S —^[F v {t)} 2 



(18) 



r p 4++;++/ 
rie /i OPI 



M H^ + Re^+ + + ^). 



F 2 4 g 



(19) 
(20) 
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The term in g _1 is due to the Coulomb photon exchanged between the scattered particles in 
diagram (d) of Fig. |I[ It should be absorbed in the static characteristics of pionium || |27|, p8| 
within the treatment of the bound state properties. Furthermore, the infrared divergent terms 
show up only in the O(q) terms. The regular terms, Re A^' and Re A^' ++ , constitute the 
central object of the current work. Their isospin limits are nothing else than 2a° + Oq and a 2 , 
respectively, 

' = i(2a° + a 2 ) + Aa (+-;+-), (21) 



32tt tnr - 6 
-LReA++; + + = a 2 + Aao(++;++) . ( 22) 

The NLO expressions for the electromagnetic corrections are found to be 

+-) = + ^ (2 + Is) - ^ (24 - , (23) 

7T 7T 7T 

9 A M 2 A _ p 2 M 2 

32 7 rAa Q ( ++;++ ) = ^ + ^| (3 + 2l 3 + 8I 4 ) - ^ (20 - , (24) 

with the combinations of the electromagnetic low-energy constants defined in (14]) and (|i~5|). 



4. For the numerical estimates of Eqs. (|23| ) and (p4|) , we use the following values |p9 
The charged and neutral pions masses are M^i = 139.570 MeV and M^o = 134.976 MeV, 
respectively. The pion decay constant is taken as F n = 92.4 MeV. The values for the strong 
sector low-energy constants are: I3 = 2.9 ± 2.4 |[23|| , Z 4 = 4.4 ± 0.3 ||3()|| . As for the two-flavour 



low-energy constants fcj, no direct numerical estimates exist up to now. On the other hand, 



some of the corresponding three-flavour counterterms Kl(p) = K\ were evaluated in |3T], |32], |33]j . 
In order to give a numerical estimate of the combinations ([14]) and (|i~5|), it is then necessary 
to relate the k\\ to the K^s (note that the combinations of fc['s involved in the amplitude 
were expressed in terms of the K^s in [12] and 16[). In the present case, this 



7T + 7T 



7r 7T 



requires to compute the corresponding three-flavour representations of the quantities Re A tb ^' 
and Re A^' |34|], and then to match them with Eqs. ( pT|) and (|22|), respectively, in the limit 
of a very massive strange quark 
between the Zf's and the L^'s, we obtain 



With the help of the strong sector matching relations 35 



K 



++;- 



[-12(K[ + K r 2 ) + 9(6*5 + Kl) - 10(^ + K$ + 72(K£ + K[ + K r u )} 



6Z 1 + In 



B m 6 
V 2 



1QZ / 4 + 21n 



M 2 



/// 



(25) 



2Z 



[12{K[ - 5K r 2 ) - 9(6K r 3 + 5Kfi + 2(5K 5 r - 31/Q + 12[K r % + K r w + K r n )\ 



In 



B Q m t 



16Z ( Z 4 — In 



12 In- 



M 2 



(26) 



where Z is the three-flavour analogue of Z, 



Using the three-flavour to two-flavour matching relations [ |35|j 

8M 2 



F 
B 



B 



Ml ln M| 



32tt 2 F 2 

M 2 
9Qn 2 F 2 



+ 



Ml 



J- n 



(2L; - L\ 



with 



Ml 



B m s 



Ml 



4 



B m s , B 



(uu) 
~F~ T 



m u =m t i=ms=0 



together with the three-flavour and two-flavour |l£| calculations of we get 



1 



32M| 
— f^ L i 



+ 



r 



Note that the following replacements are valid at the order we are working 

B m s ^Ml ± - l -Ml ± , 



with M K ± = 493.677 MeV. Using the values of the KT'a derived in § at ^ = M p = 770 MeV 
and assigning to each of them an uncertainty of ±l/(167r 2 ) coming from naive dimensional 
analysis, we obtain for the combinations (E5I) and 



e 2 M 2 . 
7r^_ tr+- 

F 2 



e 2 M\ ^ 

8.63 ±12.02, -y^lC- 



-15.13 ± 14.62. 



(27) 



This entails (we first show separately the contributions of each of the three terms in Eqs. (p3|) 
and (P|)) 



Aa (+-;+-) 
Aa (++;++) 



2.9 • 10" 3 + (1.9 ± 0.9) • 10~ 4 + (2.5 ± 7.6) • 10~ 4 
(3.2 ±0.8) • 10~ 3 , 

2.9 • 10" 3 + (8.7 ± 1.1) ■ 10~ 4 + (-12.0 ± 9.2) • 10" 
(2.5 ±0.9) • 10~ 3 . 



(28) 



From the last result, we see that the size of electromagnetic corrections to the combination 
2a[j ± Oq represents ~ 5% of the NLO value. This amounts to one-half of the size of the 
next-to-next-to-leading order (NNLO) strong interaction corrections f9j. Also, the bulk of the 
corrections arises from the electromagnetic pion mass difference, whereas the uncertainties are 
dominated by the error bars on the determinations of /C + ~ ;+ ~ and /C ++;++ . 



5. In the present work, the scattering amplitude for the process 7r + 7r~ — > n + n~ was calculated 
at LO and NLO in the chiral counting and in the presence of electromagnetic corrections. We 
have also evaluated the influence of the latter on the corresponding combination 2a(j + a 2 , of 5"- 
wave 7T7T scattering lengths which is relevant for the 2S — 2P strong energy level shift of pionium. 
The LO electromagnetic corrections come entirely from the difference between the charged and 
neutral pions masses. As for the NLO electromagnetic corrections, their size amounts to ~ 10% 
of the one for the lowest order isospin breaking effects. This corresponds to a 5% correction to 
the strong 2S - 2P level shift at NLO. 
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